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Abstract 

The recent availability of new data sets during the recent extreme solar minimum 
provides an opportunity for testing the performance of the International Reference 
Ionosphere in historically undersampled regions. This study will present averages 
and variability of topside ionospheric densities over Africa as a function of season, 
local time, altitude, and magnetic dip latitude as measured by the Coupled Ion- 
Neutral Dynamics Investigation (CINDI) Mission of Opportunity on the C/NOFS 
satellite. The results will be compared to the three topside model options available 
in IRI-2007. Overall, the NeQuick model is found to have the best performance, 
though during the deepest part of the solar minimum all three options significantly 
overestimate density. 
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22 1 Introduction 

23 The International Reference Ionosphere (IRI) is the internationally recog- 

24 nized model for calculating empirical ionospheric parameters such as density, 

25 composition, and temperature (Bilitza and Reinisch, 2008). IRI was founded 

26 as a joint project between the Committee on Space Research (COSPAR) and 

27 by the International Union of Radio Science (URSI). As new measurements 

28 in historically undersampled regions become available, it is important to con- 

29 tinuously expand the underlying database for an empirical model in order to 

30 improve its utility. One such area is the African sector; recently the addi- 

31 tion of more GPS and ionosonde stations has provided additional data (e.g., 

32 Akala et ah, 2010; Bolaji et ah, 2012). Additionally, it is important to look to 

33 the extremes of solar activity; during 2008 solar EUV dropped to the lowest 

34 observed values during the space age (Araujo-Pradere et al, 2011). 

35 This paper will present topside densities as measured from the C/NOFS 

36 satellite over the African sector for periods of both low and extremely low solar 

37 activity from 2008 to 2011. Variations with respect to local time, altitude, 

38 magnetic dip latitude, and season will be compared to the predicted values 

39 from the different topside models available in IRI-2007. 

40 2 IRI-2007 

41 The IRI-2007 model predicts the peak density ( N m F - 2 ) and the height 

42 ( h m F - 2 ) of the F -2 layer based on the CCIR maps for input geophysical indices 

43 Ri 12, the international sunspot number (Clette et al., 2007), and IG 12, the 
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44 “global effective sunspot number,” (Liu et al., 1983). The IG\ 2 index was was 

45 developed from a weighted average of ionosonde activity to estimate N m F 2 

46 based on the International Radio Consultative Committee (CCIR) maps of 

47 the ionosphere. However, unlike the international sunspot number, IG 12 has 

48 no lower limit, and during the recent extreme solar minimum is often nega- 

49 tive. Once the peak density is determined, several options may be invoked to 

50 determine the shape of the topside and bottomside profiles. 

51 IRI-2007 includes three options for calculating topside density: the IRI-2001 

52 model, the corrected topside model (herein referred to as IRI07-Corr), and the 

53 NeQuick model. IRI07-Corr uses a correction factor for the previous topside 

54 model as a function of altitude, dip latitude, and local time as calculated 

55 based on the topside sounder profiles from the ISIS and Alonette missions 

56 (Bilitza, 2004). The NeQnick model uses an Epstein-layer with an altitude- 

57 varying scale height (Radicclla and Leitinger, 2001; Coi'sson et al., 2006). This 

58 variable scale height reflects the transition of an ionosphere dominated by 0 + 

59 to one dominated by light ions at higher altitudes. 

60 Figure 1 shows the output for the three topside options for a sample profile 

6 1 for the December solstice of 2010. The solid lines represent profiles at local 

62 noon; the dashed lines are at local midnight. Note that the NeQuick model 

63 diverges from the original topside model at much lower altitudes than the 

64 corrected topside profile (IRI07-Corr). 

65 3 The C/NOFS Satellite 

ee The C/NOFS satellite is part of a space weather mission led by the US Air 

67 Force Research Laboratory to locate, understand, and predict equatorial iono- 

68 spheric scintillations (de La Beaujardiere et al., 2004). The C/NOFS satellite 
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69 


was launched in April 2008 into a 13° inclination orbit with perigee near 400 

70 km and apogee near 860 km. The record-low thermospheric density during 

71 2008 and 2009 resulted in a reduced satellite drag (Emmert et ah, 2010), and 

72 the perigee and apogee of C/NOFS slowly decayed to 399 km and 810 km 

73 respectively by August 2011. The elliptical orbit allows for a sampling of ion 

74 density over multiple 0 + scale heights of the topside equatorial ionosphere. 

75 Perigee precesses through all solar local times approximately every 65 days. 


76 C/NOFS is equipped with multiple instrument suites designed to study the 

77 ion and neutral populations and their effect on the propagation of communi- 

78 cation signals. Of interest to this study is the ion density (electron density for 

79 a quasi-neutral plasma), which is provided by the Retarding Potential Ana- 
so lyzer (RPA) from the Coupled Ion-Neutral Dynamics Investigation (CINDI) 

81 suite of instruments. The RPA technique is a well-established in situ method 

82 of sampling ion density, temperature, and composition (Heelis and Hanson, 

83 1998). 


84 Sample CINDI data is shown in Figure 2, along with the IRI-2007 predic- 

85 tions using the three topside options. Panel (a) shows three orbits from 4 Nov 
se 2010 (when perigee is at local noon); panel (b) shows the same for 6 Dec 2010 

87 (when perigee is at local midnight). In each case, perigee is marked by the 

88 vertical dashed blue line. Because the perigee of C/NOFS is near 400 km, the 
as various bottomside options will not be investigated in this study. By averaging 

90 both the measured and predicted values over all orbits for a sufficient number 

91 of days, the performance of IRI can be investigated as a function of local time, 

92 altitude, and geographic position. 
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4 Statistical performance of the topside options for IRI-2007 


94 For this study, the data will be broken into four seasons per year according 

95 to Table 1. Each season will consist of 91 days centered about the appropriate 

96 solstice/equinox, and the two equinoctial seasons remain separate. Summary 

97 plots from the March equinox of 2011 are not included in this study because 

98 the rapidly change solar flux makes it difficult to reconstruct the average iono- 

99 sphere over a 91-day period. Predicted densities from the three topside models 

100 are calculated every five seconds along every orbit, and the same averaging 

101 techniques applied to the data are also applied to the modeled values before 

102 comparison. Data beyond August 2011 will not be considered, as the definitive 

103 values of i?i 12 and IG \2 beyond this month have not been released at the time 

104 of this paper. While the average values of Riu and IG u are reported in Table 

105 1 for each 91-day season, the actual monthly values were used to run IRI-2007. 

loe The averaging techniques used here are similar to those applied in Klen- 

107 zing et al. (2011). In summary, a Savitsky-Golay filter (Savitsky and Golay, 

108 1964) is used on the CINDI data to remove small-scale irregularities, and 

109 only geomagnetically quiet times where K p is less than 3 for the preceding 
no 24 hours are used in the averages. For this study, we will only consider geo- 

111 graphic longitudes in the African sector (between 20° W and 60° E). The data 

112 will be averaged over local time (with a bin width of 2.5 hours), altitude (50 

113 km), and dip latitude (5 degrees). The orbit of the satellite does not allow 

114 for a complete coverage over all three dimensions, so certain slices through 
ns the three-dimensional space have been chosen to maximize the information 
lie conveyed. 

117 The data will be presented as mean values with a standard deviation. This 
ns reflects the day-to-day variability of the ionosphere, even under geomagnet- 
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ically quiet times. It is important to remember that IRI provides monthly 
averages of data. Previous studies have shown that there is a significant vari- 
ation of the peak density in the African sector, particularly near the dawn 
terminator (Akala et ah, 2010). A concurrent study by Simoes et al. (2012) 
of the Ionospheric Alfven Resonator (IAR) over Africa found that both the 
mean CINDI averages used in this study and the IRI predicted density profile 
produced a similarly large error when modeling the spectral response of an 
assumed density profile. However, there was a significant improvement in the 
IAR model when the mean density minus one standard deviation was used in- 
stead, indicating that the day-to-day variability is important when considering 
individual profiles. 

4-1 Variation with Local Time 

Figure 3 shows the performance of the three topside models in the African 
sector near 12° S dip latitude for the solstitial periods under consideration. 
This latitude was chosen because it has the most complete coverage over local 
time due to the orbital precession. For each season, the average ion density 
between 400 and 450 km is plotted as a black line, with the error bars rep- 
resenting standard deviation. The results from the three models are plotted 
as dashed lines. Note that at these altitudes, there is vey little difference be- 
tween the IRI-2001 topside model and the corrected topside model. Through 
the June solstice of 2010, all three topside models are high by more than one 
standard deviation for most of the day. For the December solstice of 2010 and 
June solstice of 2011, the IRI estimates are with one standard deviation of the 
mean density. In all cases, the NeQuick model performs slightly better than 
the other two options. 

Figure 4 shows the performance for altitudes between 760 and 810 km 
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145 (near the C/NOFS apogee) for the same six seasons. For these altitudes, the 
we three topside models produce significantly different predicted densities. For 
147 all seasons, the NeQuick model has the best performance, outperforming the 
ms IRI-2001 and IRI07-Corr models by a factor of 3 and 2, respectively. Note that 
149 there is still a significant overestimate (>2 standard deviations) during much 
iso of the daytime June solstices. 

151 Figure 5 shows the performance near perigee for the four equinoctial sea- 

152 sons. Note that not all local times for the March equinox of 2010 are available. 

153 In general, the NeQuick model has the best performance at this altitude. There 

154 are several notable exceptions to this (such as the post-sunset hours for March 

155 2010), but for these times all three models fall within one standard deviation 

156 of the the average density data. During the September equinox of 2010, the 

157 NeQuick estimates at nearly all local times fall within one standard deviation 

158 of the data; only the local times near dawn shown an overestimate. 

159 The densities for the equinoctial seasons near apogee (760 to 810 km) are 

160 shown in Figure 6. As in Figure 4, the three topside models predict significantly 

161 different densities. NeQuick performs significantly better than the other two 

162 options, but still overestimates by more than one standard deviation during 

163 the day in 2009. 

164 4.2 Variation with Dip Latitude 

165 The variation of the afternoon ion density with respect to magnetic dip 
lee latitude is shown in Figure 7. Only data near perigee (400 to 450 km) will be 
i67 used for the dip latitude plots in order to have the best opportunity to observe 
lee the southern Equatorial Ionization Anomaly (EIA) crest in the CINDI data. In 
i69 all cases, NeQuick provides a slight advantage in calculating the density. There 
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is little difference between the IRI-2001 and IRI07-Corr models for this region. 

171 For December 2008 (Figure 7a), the EIA crest appears in the data about 4° of 

172 magnetic latitude north of its predicted location by IRI. By December 2010 

173 (7e), the EIA crest appears in the data at or near its predicted location. 

174 The anomaly crest near midnight is shown in Figure 8 for the same six 

175 seasons. Note that for nighttime, there is a greater divergence between the 

176 IRI-2001 and the corrected topside models for dip latitudes greater than 6° S. 

177 NeQuick provides the best estimate of the three options; it performs particu- 

178 larly well for the southern half of the two equinoctial plots. For June 2010, the 

179 EIA crest appears about 3 degrees north of its predicted location at midnight. 

iso Figure 9 shows the afternoon position of the southern EIA crest for the 

181 periods of extremely low solar activity. Note that the position of the anomaly 

182 crest is unclear in all but September 2010 (9d). The nighttime plots with 

183 respect to dip latitude are shown in Figure 10. 

184 5 Discussion 

185 For all seasons and solar activities, the NeQuick topside option performs 
lee best when compared to average in situ measurements of ion density. This is 
i87 consistent with previous studies of vTEC over Africa (Adewale et ah, 2012), 
we as well as topside sounder data from the ISIS and Alouette missions (Bilitza, 

189 2009). Additionally, the NeQuick model represents a significant improvement 

190 over the other options at high altitudes. Indeed, the performance at the higher 

191 altitudes is likely to be the cause of the discrepancy previously reported in 

192 vTEC measurements (Adewale et ah, 2012). Note that there are some seasons 

193 where NeQuick is more likely to overestimate at higher altitudes (in particular, 

194 the June solstice), indicating that the model does not decrease with altitude 



195 


as rapidly as it should. 


we There tends to be an overestimate near the dawnside terminator between 4 

197 and 8 SLT for altitudes near the F-peak, which goes away for higher altitudes. 

198 This is consistent with observations from GPS vTEC over Africa (Adewale et 

199 al., 2012) and observations from CHAMP in the equatorial region (Lee et al., 

200 2011 ). 

201 During the recent solar minimum between cycles 23 and 24, solar EUV flux 

202 dropped to record low levels, well below what would have been expected based 

203 on observations of F10.7 (Chen et ah, 2011). This was confirmed not only from 

204 direct measurements of EUV (Araujo-Pradere et al, 2011), but also through 

205 the demonstration that the neutral atmosphere was contracted beyond expec- 
2 oe tation through satellite drag calculations (Solomon et ah, 2010) and direct 

207 measurement of the neutral scale height (Haaser et ah, 2010). Ionospheric be- 

208 havior during the extreme solar minimum was more complicated than that 

209 of the neutrals because the response to EUV changes is non-linear (Araujo- 

210 Pradere et al, 2011). The vertical ExB drifts were found to be significantly 

211 different, including the absence of the pre-reversal enhancement during 2008 

212 and 2009 and downward drifts during the early afternoon (Pfaff et ah, 2010; 

213 Stoneback et ah, 2011). 

214 Previous studies utilizing the CHAMP and GRACE satellites have shown 

215 that IRI-2007 topside overestimation occurs as early as 2006 (Liihr et ah, 

216 2010; Lee et al., 2011); this study along with the global equatorial averages 

217 presented in Klenzing et ah (2011) shows that as solar cycle 24 proceeds, the 

218 performance of the IRI model — specifically, NeQuick — improves significantly. 

219 The El A density crest, sometimes referred to as the Appleton anomaly, 

220 is a result of the equatorial fountain effect, where plasma near the magnetic 
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equator is lifted up by ExB drifts and slowly diffuses down the field lines 
(e.g., Schunk and Nagy, 2000). During December of 2008, the EIA crest as 
observed in the CINDI data appears about four degrees closer to the magnetic 
equator than predicted. This discrepancy is likely to be an effect of the con- 
tracted ionosphere, and may be the result of the anomaly crest moving down 
in altitude, moving in to neighboring flux tubes, or a combination of both. 


6 Summary 


The performance of the three topside model options of the International 
Reference Ionosphere (IRI) were compared to in situ measurements of density 
from the C/NOFS satellite in the African sector (20° W to 60° E) for periods 
of low and extremely low solar activities. The major hirelings were: 

(1) The NeQuick topside model performs the best overall when compared 
to the mean in situ data. The only exceptions to this observation occur 
when all three models are within one standard deviation of the mean. 

(2) At higher altitudes (~800 km), there is a significant overestimate of den- 
sity by the IRI-2001 and IRI07-Corr topside models for all seasons com- 
pared. 

(3) During the extreme solar minimum of 2008 and 2009, all topside options 
tend to overestimate density in the topside ionosphere over Africa. 

(4) The southern EIA crest over Africa is about 4 degrees closer to the mag- 
netic dip equator than predicted by IRI during the December solstice of 
2008. 

(5) The performance of IRI in the summer hemisphere is dramatically im- 
proved during 2010. IRI still tends to overestimate electron density during 
the winter for dip latitudes greater than 15 degrees. 
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Table 1 


Seasonal divisions for the topside profile reconstructions, including the average solar 


activity represented by F10.7, Rii 2 , and IG 12 for each period. 


Period 

Range of Days 

F10.7 

CtFlO.7 

R 112 

IG 12 

Dec Solstice 2008 

6 Nov 2008 

- 4 Feb 2009 

67.2 

0.9 

2.0 

-9.2 

Mar Equinox 2009 

3 Feb 2009 

- 4 May 2009 

69.1 

1.2 

2.0 

-8.4 

Jun Solstice 2009 

7 May 2009 

- 5 Aug 2009 

70.1 

2.0 

3.2 

-6.3 

Sep Equinox 2009 

8 Aug 2009 

- 6 Nov 2009 

70.7 

2.5 

6.2 

-0.8 

Dec Solstice 2009 

6 Nov 2009 

- 4 Feb 2010 

75.1 

4.7 

8.5 

3.9 

Mar Equinox 2010 

3 Feb 2010 

- 4 May 2010 

80.8 

5.0 

12.6 

10.4 

Jun Solstice 2010 

7 May 2010 

- 5 Aug 2010 

77.7 

5.4 

16.8 

15.7 

Sep Equinox 2010 

8 Aug 2010 

- 6 Nov 2010 

81.2 

4.4 

21.1 

19.1 

Dec Solstice 2010 

6 Nov 2010 

- 4 Feb 2011 

81.1 

4.2 

29.4 

26.2 

Jun Solstice 2011 

7 May 2011 

- 5 Aug 2011 

98.8 

10.6 

53.4 

56.1 
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Fig. 1. Sample IRI-2007 profiles, illustrating the differences between the three top- 
side models. Density profiles at local noon are solid lines; profiles at local midnight 
are dashed. 
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Fig. 2. Sample CINDI data for three consecutive orbits, compared with the IRI pre- 
dicted densities for each topside model. The vertical dashed lines represent perigee 
(~ 400km). 
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Fig. 3. Comparison of the three topside options of IRI for ion density with the 
CINDI data for a period of low solar activity. The data is averaged between 400 
and 450 km altitude at 12° S magnetic dip latitude. Note that the IRI-2001 and the 
IRI07-Corr models significantly overlap at these altitudes. 
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Fig. 4. The same as Figure 3, but for altitudes between 760 and 810 km, 
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Fig. 5. The same as Figure 3, but for the equinoctial seasons. 
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Fig. 6. The same as Figure 4, but for the equinoctial seasons. 
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Fig. 7. Comparison of the three topside options of IRI for ion density with the 
CINDI data for a period of low solar activity. The data is averaged between 400 
and 450 km altitude and between 13.5 and 16 hours of solar local time. 
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Fig. 8. The same as Figure 7, but near local midnight. Data is averaged between 
22.5 and 1 hours of solar local time. 
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Fig. 9. The same as Figure 7, but for the equinoctial seasons. 
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Fig. 10. The same as Figure 8, but for the equinoctial seasons. 
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